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We study a generalization of the Kondo model in which the impurity spin is represented by 
Abrikosov fermions in a rotation group SU(P) larger than the SU(JV) group associated to the spin 
of the conduction electrons, thereby forcing the single electronic bath to underscreen the localized 
moment. We demonstrate how to formulate a controlled large N limit preserving the property of 
underscreening, and which can be seen as a "dual" theory of the multichannel large N equations 
usually associated to overscreening. Due to the anomalous scattering on the uncompensated degrees 
of freedom, the Fermi liquid description of the electronic fluid is invalidated, with the logarithmic 
singularities known to occur in the S = 1 SU(2) Kondo impurity model being replaced by continuous 
power laws at N = oo. The present technique can be extended to tackle the related underscreened 
Kondo lattice model in the large JV limit. We discover the occurence of an insulating pseudogap 
regime in place of the expected renormalized metallic phase of the fully screened case, preventing 
the establishement of coherence over the lattice. This work and the recent observation of a similar 
weakly insulating behavior on transport in CeCuAs2 should give momentum for further studies of 
underscreened impurity models on the lattice. 



I. INTRODUCTION 

The Kondo model, initially introduced to describe with 
great success the behavior of diluted magnetic impurities 
in metals Q], has been an exciting playground for con- 
densed matter physicists in the last two decades. Indeed, 
the extreme simplicity of this many-body Hamiltonian 
was key to some unexpected progresses in the field, such 
as the discovery of the effective theory describing the 
Mott transition @, and as a prototypal description of 
Non Fermi liquid behavior in multi-channel [j| and pseu- 
dogap extensions Localized quantum impurities also 
provide the building entity to describe heavy fermion ma- 
terials Q , and this is still a strong motivation for study- 
ing Kondo models on the lattice, most interestingly in 
relation to the violation of Fermi liquid behavior asso- 
ciated to the breakdown of Kondo screening close to a 
quantum critical point [(|, as discovered in a wide class 
of /-electron metals such as CeCue-xAiia 0. Interest 
in the Kondo model was also revived by the possibility 
of building and controlling precisely quantum impurity 
states in semiconducting devices, leading to the realiza- 
tion of the Kondo effect in quantum dots Q • The possi- 
bility of observing Non Fermi liquid behavior in a multi- 
channel setup was recently advocated 0, 0, and 
provides further incentive for a complete understanding 
of quantum impurity problems. 

Quite recently, peculiar attention was give n to the 
Kondo problem by Coleman and Pepin [l^. in the 
study of an underscreened model which showed relatively 
strong deviations from Fermi liquid behavior in the spe- 
cific heat (we point out however that related effects can 
be traced back to the physics of the ferromagnetic Kondo 
model 01 1 see also This surprising result was inter- 

preted as due to anomalous scattering of conduction elec- 
trons onto the remaining unscreened spin degrees of free- 



dom, giving rise to a "singular" Fermi liquid fixed point, 
which differs from the intermediate fixed point 0, 0] as- 
sociated to usual Non Fermi liquid by the fact that the 
effective ferromagnetic Kondo coupling renormalizes to 
zero at low energy Q. This insightful work was how- 
ever limited to a range of frequencies higher than the 
Zeeman energy, and concentrated mainly on thermody- 
namic quantities. In order to study in more detail the 
dynamics due to these anomalous excitations in the un- 
derscreened Kondo model, we would like to find a simple 
approach that would be able to grasp the full crossover 
from the local moment regime at high temperature down 
to the underscreened state in which nearby conduction 
electrons are tightly bound to the impurity. In the real- 
istic case of SU(2) spins, we note that the deviation from 
exact screening can be in principle tuned Q by chang- 
ing the size 5 of the impurity spin and/or the number of 
screening channels M, obtaining a transition from under- 
screening at 25 > M to overscreening at 25 < M (with 
perfect screening at 25 = M). 

A simple route to capture this physics consists in es- 
tablishing a large N limit of the problem, generalizing 
the quantum spin to the SU(iV) group. A crucial step 
lies however in the choice of the SU(iV) representation in 
which the spin is considered. It is known that fermionic 
representations at large N only allow perfect scree ning 
when M <^ N ^(| and overscreening otherwise [rH Il8j . 
Interestingly, bosonic representations of SU(iV) preserve 
the distinction between "small" and "large" spin as found 
by Parcollet and Georges [l^ > and ahow to study the un- 
derscreened situation. However, this case presents some 
pathologies: the T-matrix scales as 1/N, and the singu- 
lar behavior as found in [l^ seems to be absent of the 
solution poj . 

The previous discussion illustrates the need for an al- 
ternative large N limit to describe the underscreened 
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Kondo effect, but also gives momentum to the idea we 
will pursue in the following. Indeed, we can under- 
stand that overscreening and underscreening are some- 
what "dual" in the sense that, while the former situa- 
tion is reached in the presence of many screening chan- 
nels, the latter case should be obtained by considering 
many "spin channels" . A simple way to formalize this 
is to strongly enlarge the symmetry group of the impu- 
rity spin, thus forcing underscreening by the bath of con- 
duction electrons. In fact, if one considers a generalized 
Kondo model involving a single bath of SU(TV) electrons 
interacting with a localized SU(P) spin, where P is larger 
than TV, one expects underscreening, independently of the 
representation chosen for the impurity. Because the M- 
channel large TV limit which describes overscreening is 
obtained for M oc N |l7l [l8| , we can guess that a reason- 
able large TV limit of underscreening can be found when 
P = KN (K oc TV being the number of "spin channels" ) . 
We will show that the saddle-point equations derived in 
this way at TV = oo present a "dual" structure to the Non 
Crossing Approximation associated to the overscreened 
case 0[ 0, , and reveal non Fermi liquid behavior 
with a continuously varying exponent parametrized by 
the ratio 7 = P/N 2 . The lattice version of this gen- 
eralized underscreened Kondo model is also straightfor- 
wardly solvable, due to the absence of generated RKKY 
interactions. Because singular scattering of delocalized 
electrons on the degenerate manifold of unscreened spin 
occurs, coherence cannot establish over the lattice, in 
contrast to the exactly screened Kondo lattice [12 ■ This 
proves the existence of a generic insulating- like state with 
a pseudogap density of states, parametrized by the same 
anomalous exponent unveiled in the single impurity un- 
derscreened Kondo model. We conclude the paper by 
future prospects regarding the theory of underscreened 
Kondo models, and by the possible observation of their 
physical consequences in strongly correlated materials, 
pointing out some similarities with the cerium based com- 
pound CeCuAs2 , on which transport measurements were 
recently reported. 

The remainder of the paper is organized as follows: 
in section [H]wc demonstrate how the new large TV limit 
of the underscreened Kondo impurity model can be per- 
formed, followed by a physical discussion of the results 
in section ITTT1 The Kondo lattice extension is examined 
in section HV1 



1. Heuristic derivation 

The basic idea, which was presented in the introduc- 
tion, is to consider a single-channel Kondo Hamiltonian 
involving conduction electrons with TV spin flavors and 
interacting with a SU(P) spin S mm i localized at the ori- 
gin: 

ff = 5>4/W+ £ J™P'c\ a c kla ,S m , m (1) 

ka kk' crcr' mm' 

where a = 1 . . . N, m = 1 . . . P. The matrix of coupling 
constants J™, n will be specified in the following. We 
choose a completely antisymmetric representation of the 
spin using Abrikosov fcrmions: 

S m 'm — fm'fm ~ ~p^mm' (2) 

Q = J2f™fm (3) 

m 

where the second equation is the necessary constraint to 
enforce the spin size Q = q$P (this scaling, where go is of 
order 1, is necessary to get a large TV limit, see below). As 
emphasized previously, P > TV leads to underscreening, 
and a likely scaling to obtain a solvable limit is P oc TV 2 . 
Let us therefore set P = KN, with K = 7 TV, expressing 
fcrmions in a double index notation: 

S m 'm — fa'a'faa ~ "J^^aa'^aa' (4) 

where m represents the doublet of indices (a, a), with 
a = 1 ... K. Neglecting potential scattering terms, which 
only contribute to next order in 1/TV, the Hamiltonian 
now reads: 

h = J2^ c I c ^+J2 (5) 

ka kk' aa' aa' 

This is now completely general, and we would like to 
chose the simplest coupling between the itinerant SU(TV) 
fcrmions and the localized SU (KN) spin. In the spirit of 
having K "spin channels" , we set: 

jaa'aa' _ f fi \ 

J oV — AT aa V ) 



II. NOVEL LARGE N LIMIT OF THE 
UNDERSCREENED KONDO MODEL 

A. Model 

Our aim in this section is to introduce a Kondo model 
in which the underscreened aspect is built in from the 
beginning, while a simple large TV limit of the problem 
can be found by proper rescaling of the parameters. 



and obtain the Hamiltonian 

H = Ys ekC kcr C k<T+J2 ^ C L C k'a'fL'faa (7) 
ka kk'aa'a. 

loP = £/JU«r (8) 

aa 

associated with the single constraint eq. © ■ This will be 
our starting point for the large TV solution. 
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2. Alternative interpretation 

One might argue that coupling a band of fermions to 
a spin with a different symmetry group is not a very 
transparent concept on a physical point of view. Indeed, 
one can try to view alternatively the Hamiltonian (J7J) 
as the standard SU(iV) Kondo model with a localized 
SU(7V) spin taken in a rectangular representation 23] of 
size K x qoN: 



5' 



/ j J acrJ ac 



(9) 



Usual Kondo coupling to the SU(iV) electronic bath leads 
indeed to the Hamiltonian Q ■ We will however not pur- 
sue this route in the present work, as this type of repre- 
sentation asks to consider a set of N 2 constraints: 



qoN5 aa < = fLfa'c 



(10) 



which makes the large N limit impracticable. We will 
come back to the differences between the large group 
technique (developed in this work) and the rectangu- 
lar representation approach in section fill C 21 while dis- 
cussing the physical results. 



B. Derivation of the saddle-point equations 

In order to solve the model ^ 718(1 . we first write the 
imaginary time action of the problem with inverse tem- 
perature /3: 

s = A[e/-(^ + a )/-+E c L(^ + ^) c 

J ° 1 » ka 
r@ r j 

+ / dT M H C L c fc' CT ' /<Lr'/aa-9oPA 

•* L 7,7../ /~ 



kk' aa' a 



introducing a Lagrange multiplier A to enforce the con- 
straint . The precise form of the action is hinting 
that a large N solution is possible, which we perform 
now. 

The following step is to decouple the interaction with 
a bosonic field : 

L a kaa v 



and integrate out the fermions y aa : 

S = [dr\j2cUdT + e k )c k „+J2^T l -<loP* 
Jo I , J 



k 

/dr /dr'l G/ o(r-r') £ (c^Bj^flt^)^ 



(13) 



where Gfo(iuj n ) = l/(iu> n — A), with u> n — (2n + l)ir/ (3. 
We finally introduce the electron sitting at the impurity 
center, c\ =J2k C L> so tnat: 



S = /dr £^-, PA 
Jo I a J 

fdrfdr'G-^r-r^cUryAr') 
Jo Jo a 

+ f%r [c\r'^G fQ (r -r')J2(ciB a ) T (Bi C<7 ) T , 
Jo Jo ^ 



(14) 



where G c o(iuj n ) = \/{iu> n — e k ). Because K = jN 
scales as N, the existence of a saddle-point is manifest 
in the previous expression. Following |18|. we find the 
integral equations: 



G c {iuj n ) = (cl(iu! n )c a {iuj n )) = — — j- 



1 



G c0 (iuj n ) - Y, c {iuj n ) 
(15) 

G B {iv n ) ee (Bl{iv n )B a (iv n )) = 1 (16) 

S c (r) = 7 G f0 (T)G B (T) (17) 
E b (t) = G f o(0-T)G c (T) (18) 
9o = - 0+)/ QCT (r = 0)) (19) 

where v n = 2mr/f3 is a bosonic Matsubara frequency. 

C. Interpretation of the formalism 

Let us first comment on the formal analogy that the 
system of equations (I15I18|) shares with the usual Non 
Crossing Approximation (NCA) [13, EE EH- We notice 
indeed that the roles of the local bath electron cj. and 
of the Abrikosov fermion f^ a are simply exchanged with 
respect to the usual NCA structure, which, we recall, is 
associated to overscreening ^3- In particular, the self- 
energy now involves the propagator Gfo(iuj n ) — l/(iu) n — 
A) instead of G c o{iuj n ) — J2k V(* w n — £ fe) m the NCA, 
and this difference will be shown below to radically affect 
the physics, with the occurence of underscreening instead 
of overscreening. This interesting "duality" is one of the 
main results of the present work. 

It is also very appealing to remark that the computa- 
tion of physical observables related to the itinerant band 
does not involve directly pseudo-fermion degrees of free- 
dom, contrarily to the multichannel case. Indeed, the 
present theory is formulated directly in terms of the bath 
propagator, to which the T-matrix is simply related by 
the relation: 

G c (k,k',iuj n ) EE (c\ l(r {ioJ n )c ha {iu) n )) (20) 

= V T(iu n )- 

iuj n - e k VjJ n - e k iu) n - e k > 

(21) 
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If we sum the previous expression over all momenta, we 
find: 



G c (iuj n ) — G c0 (iuj n ) 

G 2 cQ (iuj n ) 



1/E c (iuj n ) - G c0 {iuj n ) 



(22) 
(23) 



where (|15|) was used to obtain the second expression. In- 
terestingly, we note that the T-matrix l|23|) is of order N° 
in the present scheme, whereas it happens to scale as 1/N 
in the multi-channel large N limit a certain draw- 
back for the theory of overscreening. However, by the 
"duality" argument, we expect to find that the spinon 
propagator Gf(iw n ) = (fl a (iuj n )f a<J (iu! n )) only shows 
1/7V corrections to the free impurity limit (J = 0), which 
can indeed be easily checked by a direct computation: 




FIG. 1: (color online). Logarithmic plot of 
7 = 1,3, 10 (bottom to top) at /3 = 1000, J ; 



Gf(u> n ) — G /o(to n ) — G j {iuin)G x (iuin) (24) lines are fits to a power law form |w| 



TrriE c (iUn) for 
1, D = 4. Thin 
according to Eq. I31H . 



G x (r) 



-G c (t)G b (-t). 



(25) 



This means that, although the T-matrix is conveniently 
captured by the present scheme, computation of e.g. the 
spin susceptibility involves necessarily 1/N corrections. 
This remark allows however to solve explicitely the con- 
straint equation l|19|) at the leading order: 



qo = G /0 (r = Q-)+O(i) 
= ilogi^ ^ + 



l 



S /3A + i 



(26) 
(27) 



where C c and Cj, are undetermined constants. This ana- 
lytical result is well borne out by the numerical solution 
of the saddle point equations as shown by figuren(in the 
following computations, a semi-circular density of states 
with half- width D was chosen to model the bath of con- 
duction electrons). 

We would like to contrast this peculiar low energy be- 
havior with the one observed in a standard Fermi liquid 
(as in the fully screened Kondo model). The interpreta- 
tion of these results can be made clearer by considering 
the T-matrix of a related Anderson model: 



III. PHYSICAL STUDY OF THE NON FERMI 
LIQUID REGIME 

A. At particle-hole symmetry 

For simplicity, we start with the assumption of particle- 
hole symmetry, qo — 1/2 (the free bath G c o(io; Tl ) will 
always be assumed symmetric in the following). This 
implies that the constraint l|19|) is fulfilled provided A — 
0, and we have: 



G/o(r) 



1 



Sgn(r) 



i.e. long range correlations in the self-energies I|17ll8fl . 
By analogy with the low-frequency analysis usually per- 
formed in studying the NC A equations 0, 0] , we as- 
sume power law behavior of the self-energies and plug 
into the saddle-point equations. Self-consistency can be 
achieved at zero temperature (see Appendix and we 
obtain for real frequency quantities: 



Tm Sb(w) 



C B M+ Q SgnH 

2 1 
— arctan 

7T y/j 



(29) 
(30) 

(31) 



T(iu) n ) 



V 2 



iLOn - T, d (itUn) - V 2 G c o{iuj n ) 



(32) 



where the local fermion dj. has an hybridization V to 
the cJ.-electrons and is subject to a local Coulomb re- 
pulsion term U (and possibly to a strong Hund's rule if 
one wants to stabilize aS = l spin in order to make the 
connection to the underscreened case), giving rise to the 
term Y,d{iu n ) in the previous equation. Comparing ll'.'il) 
and H32|) . we see that the inverse self-energy of the bath 
electrons is simply related to the impurity self-energy by 
1/Y, c {bj) = [u>— T,d{ijj)]/V 2 . In a Fermi liquid, we know [l| 



(28) that the impurity self-energy obeys at low frequency: 



Sf L (w) = (1 - 1/Z)u> + iAuj 2 + 



(33) 



where Z is the quasiparticle residue and Au> 2 the inelastic 
scattering rate. The previous general identification be- 
tween 1/E C and provides the self-energy of the local 
electron in the bath for a Fermi liquid: 



Y 



ZV 2 



Sf(a 



+ iAV 2 Z 2 



(34) 



the first term corresponding to regular elastic scattering 
from the impurity. The fact that this self-energy diverges 
as a power law at low frequency in the underscreened 
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FIG. 2: (color online). Imaginary part of the T-matrix for 
decreasing temperatures /3 = 10, 100, 1000, oo (bottom to top) 
with J = 1, D = 4, 7 = 1. The left inset is a zoom on the 
Kondo resonance, and the right inset shows the associated 
depletion of low energy states in the local spectral function of 
the bath (-l/ir)lmG c (u>) oc \uj\ a . 



particle-hole symmetry is broken. The propagator ap- 
pearing in the bubble giving the self-energies, equa- 
tions I(17I18|) . is given by: 



G f0 (r) = — 



A/3 + 1 



for < t < (3 



(37) 



and decays exponentially. We argue that the long range 
correlations which are crucial to maintain the non-trivial 
power laws are still present, because f3\ saturates at low 
temperature, from 1(27(1 . This would imply that the non 
Fermi liquid state survives the introduction of the asym- 
metry parameter go 7^ 1/2, as can also be verified from 
the numerics. 

However, the structure of the saddle point is such that 
the constraint © scales as A 2 instead of N, and hints 
that it might be necessary to account for 1 /N corrections 
to the result 127(1 . We do not find that this possibility 
really modifies the previous result, A oc 1//3, although we 
think this question deserves further clarification. 



case 1(29(1 signals anomalous scattering on the remaining 
unscreened degrees of freedom, which ultimately violates 
the Fermi liquid description of the problem. This result 
can be more clearly witnessed from the impurity self- 
energy itself 



r 2 



E c (w) 



ex M Q 



(35) 



showing a power law scattering rate for the localized state 
which even dominates the linear part found in the self- 
energy of a Fermi liquid, equation J33J). Indeed, we de- 
duce from this behavior a frequency dependent quasipar- 
ticle residue Z(uS) = [1 — <9£d(o;)/<9u;] -1 cx M 1-a which 
vanishes at low energy, demonstrating a clear breakdown 
of Fermi liquid theory. 

A final way to put into evidence the singular energy 
dependence of this underscreened impurity is to compute 
the T-matrix. We find indeed that the Kondo resonance, 
whose development is illustrated on figure [3 displays a 
cusp at low frequency: 



-XtoT(w) 



7T/O0 



flu 



(36) 



using in equation 1(23(1 the low frequency behavior of the 
free Green's function, G c o{uj) — —iirpo. This expression 
shows that, although the unitary limit is recovered in the 
T-matrix at w = 0, Fermi liquid behavior is nevertheless 
violated by the non-analytic expansion a low frequency. 



B. Away from particle- hole symmetry 

Here we discuss briefly the case q ^ 1/2, and inves- 
tigate whether the previous results remain valid when 



C. Discussion: theoretical aspects 

1. Possible extensions of the formalism 

The present approach seems well suited to describe 
the underscreened situation of quantum impurity mod- 
els. We can however point out two kinds of limitations, 
which are in fact also inherent to the Non Crossing Ap- 
proximation of multichannel models as well |18j . First, 
although the T-matrix can now be naturally extracted 
at the saddle-point level, some physical quantities such 
as the spin susceptibility appear only at next to lead- 
ing order in X/N. Second, motivated by some puzzles 
raised in heavy fermion compounds, it remains of great 
interest to develop faithfull large A expansion methods 
that can work in the exactly screened case too (an inter- 
esting step in this direction was however taken in [19|). 
It was in fact already acknowledged [Hj that a simple 
generalization of the Kondo model is quite suitable for 
overcoming both of these difficulties: it consists in tak- 
ing the localized spin in a rectangular representation of 
SU(iV), as indeed done in the present work (if one were 
able to take good care of the necessary constraints, see 
discussion in section III A 2(1 . while assuming also a large 
number (of order A) of screening channels (in this paper 
we only considered the single channel case). This ex- 
tension of the formalism, also known as "matrix Kondo 
model" , leads however to technical difficulties due to pro- 
liferation of Feynman graphs at large A. It is actually 
possible to resum this complicated diagrammatics u sing 
a soft constraint approach inspired by a recent study [lOj 
of the interplay of Kondo effect and Coulomb blockade 
(viewed now as a purely mathematical trick). Such a 
theory would be a bridge between the usual large A ap- 
proach of overscreening and the present work, both ap- 
pearing as simple limitating cases of it, thus finding an 
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unity behind our "duality" picture. Such a theory will 
be considered in a forthcoming publication [2(j 



IV. PSEUDOGAP FORMATION IN THE 
UNDERSCREENED KONDO LATTICE 



A. Model and large TV solution 



2. On the nature of the low energy singularities 



As we have demonstrated from the analytical and nu- 
merical solution of the new large-N equations, the un- 
derscreened Kondo impurity model displays violations 
of low energy Fermi liquid behavior, although those are 
less severe as the one observed in the overscreened case 
(where the deviation from the unitary limit associated to 
finite lifetime of the electrons at zero temperature even 
dominates the anomalous frequency dependence of the 
scattering rate). The underscreened impurity is never- 
theless associated to Non Fermi liquid behavior: the lo- 
cal quasiparticle weight vanishes at low energy and the 
system presents additional degeneracy due to the non 
fully screened spin. It is interesting to compare this re- 
sult to previous exact knowledge from the Bethe ansatz 
solution of underscreened Kondo models. Although the 
general features we have sketched so far are indeed ob- 
tained from the Bethe ansatz analysis, one crucial differ- 
ence concerns the precise nature of the singular frequency 
dependence of the physical quantities: while we observe 
universal non-trivial exponents in the large N solution, 
the 5 = 1 SU(2) Kondo model rather displays logarith- 
mic behavior [12j (although power laws are actually found 
in the anisotropic case [lj])- We think this difference is 
a minor point, that we would naively like to attribute to 
special features associated to the large N generalization 
of the model. This point is actually more subtle. In- 
deed , the hamiltonian we have solved can be argued 
(see section III A 2|) to be equivalent to a SU(iV) Kondo 
model with spin in a rectangular representation, which is 
actually diagonalized by Bethe ansatz [27|. and found to 
display also logarithmic singularities! As stressed previ- 
ously, we have however enforced a different type of con- 
straint in the present approach, and we would ultimately 
pin point this difference to be the origin of the discrep- 
ancy in the low energy behavior as compared to the exact 
solution. This discussion would tend to illustrate in any 
case the fact that the physics can be sensitive in how 
the constraint is implemented, but we will reinforce our 
judgment that the present method captures the essential 
features of underscreening. As we will now demonstrate, 
the main interest of our approach is that it can apply to 
problems where no exact solution is available, in particu- 
lar in the many impurities extension of the underscreened 
Kondo model. 



We introduce here the lattice extension of the previous 
single impurity model, which consists of a dense network 
of impurities carrying a SU(P) spin on which a band of 
SU(iV) itinerant electrons scatter: 



- V ct 

jcra' a 



C ja' fjaa'fjaa (38) 



Here j labels sites, with 



-kcj 



= 4 e ' M] > v is tlic 



chemical potential in the c-band and other conventions 
are similar to previously. This model is not very real- 
istic in many respects, and to date lacks a physical re- 
alization in condensed matter systems (see however the 
ending discussion in section liVDI) . On a purely theoret- 
ical point of view, we can raise however two interesting 
questions. First, what is the behavior of the electronic 
band below the Kondo temperature associated to under- 
screening? For this to make sense, we need to argue 
that the magnetic processes involving the remanent spin 
degrees can be ignored, for example due to frustration 
(maybe by analogy to the compound L1V2O4). Second, 
when the partially-screened local moments start to or- 
der, and quench their macroscopic entropy, how does the 
system behaves? In the following large TV analysis of this 
problem, we will see that intersite magnetic correlations 
are absent, and that the first of these questions can be 
answered. 

The large N limit is derived following the same steps 
performed in section [H] introducing local bosons B^ a to 
decouple the Kondo term on each site, and integrating 
the Abrikosov fermions: 



Jo , 

u krr jq 



k 

We can simply read off this expression the final saddle 
point equations, which are completely identical to 1)161 
I18fl . except that the local bath propagator now reads: 



J 



(39) 



G c {iuj n ) = (cl a (iu} n )c ia (iuj n )) 

= V- 

z — t 1, 



e k +/i - S c (iw n ) 



(40) 
(41) 



B. Interpretation 



The new system of integral equations I|ltjll8[l , (|41|l is re- 
markably simple in its structure, a fact due to the absence 
of intersite correlations at this level of approximation. 



7 



Indeed, as the action (|39|l shows, no magnetic RKKY in- 
teraction is generated at the saddle-point level, which is 
expected since the additional quantum number a carried 
by the localized spins cannot be transported from site 
to site by the itinerant fcrmions. We point out however 
that this is not a general feature of underscreened mod- 
els, and one can easily check that the RKKY interaction 
is indeed mediated by the electronic bath in the finite S 
SU(2) case. 

Besides, expression Q41JI is reminiscent of a self- 
consistent T-matrix approximation and signals that elec- 
trons in the bath are still anomalously scattered by the 
localized spins, acting independently of each other. In- 
deed, if we assume that S c (iw„) is divergent as in the 
single impurity case, equation Ij29(l . we see that this sin- 
gular self-energy dominates the fc-summation in l|41l) at 
low energy, so that: 



G c (iuj n ) 



1 



-S c (iw n ) 



(42) 



with the same exponent (|3 1|) as found previously. This 
means that in the underscreened Kondo lattice, there is 
no important distinction between the cases of dense and 
diluted impurities (on the point of view of the bath elec- 
trons that feel the Kondo interaction), and that the itin- 
erant electron density of states always shows a pseudogap 
at low energy. This is quite different from the situation 
of exactly screened models, where a hard hybridization 
gap would open at half-filling. In the perfectly screened 
case, coherence can always be re-established upon dop- 
ing, despite the fact that each individual impurity scat- 
ters strongly the electrons. The result l|42|) would how- 
ever let us think that electronic degrees of freedom always 
remain confined in the underscreened Kondo lattice, as 
we now check on the numerical solution of the saddle- 
point equations. 



C. Results 

We will be again interested in the T-matrix, which is 
related to the (translation invariant) c-electron Green's 
function G c (k,iu n ) = {& krT {iu>n)c krT (iu>n)) by: 



G c (k,iuj ri 



1 T(k,iu) n ) 

iui n — £fc + A* H - £fe + A*) 2 



(43) 



This relation can be best understood from an equiva- 
lent periodic Anderson model, where T{k, iLo n ) is pro- 
portional to the momentum- and frequency-dependent 
Green's function of the localized electrons. From the ef- 
fective action (|3l?)l we have: 



G c (k,iuj n ) 



1 



iu>„ - ek + (j, - H c (iu> n ) 
so that the T-matrix can be expressed as: 

1 



T(k 7 iu n ) 



1/E c (iu n ) - 1/ (iu> n - e k + //) 



(44) 



(45) 




FIG. 3: (color online). Imaginary part of the T-matrix at 
/3 = 100 and for two values of the chemical potential /i = 
(symmetric curve) and /i = 1 (asymmetric curve), with J = 1, 
D = 4, 7 = 1, showing the presence of an hybridization pseu- 
dogap (see zoom in left inset). The right inset shows the 
corresponding dip in { — l/Ti)ImG c (uo). Because of finite tem- 
perature effects, spectral weight is present at zero frequency. 



Again, we identify 1/E c (ia> n ) as the impurity self-energy 
Ed(?w n ) (up to a factor V 2 ). The local T-matrix, 
T(iu n ) = YlkT{k,iuj n ) is easily calculated from l|43|l af- 
ter the numerical solution of the saddle-point equations, 
and is shown in Figure[3](here also a semi-circular density 
of states for the c-electron was taken) . 

From the numerical solution, and in agreement with 
the previous analytical analysis, we see that an hybridiza- 
tion pseudogap opens in the spectrum, irrespective of 
the filling of the c-band. This prevents coherence to be 
reached over the lattice at zero temperature (note that 
Figure |3 is performed at finite temperature, so that the 
pseudogap in the density of states is filled by thermal ex- 
citations). Therefore, the strongly underscreened Kondo 
lattice is strictly speaking an insulator, with a resistiv- 
ity likely showing a power law increase as temperature is 
lowered, instead of the activated or metallic behavior ex- 
pected respectively in usual Kondo insulators or in heavy 
fermion metals. In order to shed some light onto this re- 
markable result, we put it into the perspective of previ- 
ous theoretical work and experimental measurements on 
a newly discovered /-electron material. 



D. Physical discussion 

1. Underscreened Kondo lattice: status of theory 

The main physical result of the paper is that the under- 
screened Kondo lattice is an insulator with a pseudogap 
density of states above the ordering temperature of the 
partially screened moments, for all filling of the electronic 
band. Theoretically, we can ask whether this behavior is 
peculiar to the large ./V solution performed here or can 
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be expected to be generic. To answer this question, it is 
interesting to notice two facts. First, previous work by 
K. Le Hur in a one dimensional model of underscreened 
Kondo ladder |28| demonstrated infact similar insulating 
behavior (this quite peculiar model was chosen in order to 
prevent the Haldane spin gap associated to the ID S = 1 
chain, which would pollute the interesting physics). Sec- 
ond, our large N solution displays actually some high 
dimensional features 0, such as a local self-energy, as 
seen in equation 144|) . Because this insulating behavior 
is found in two such extreme cases, we could expect for a 
continuity of this non-metallic behavior of underscreened 
lattice in intermediate dimension also. 

We now turn to the nature of this insulating state. As 
discussed in the single impurity case, section ITlI C 21 the 
large N solution predicts that the dynamics is governed 
by universal power law behavior rather than the loga- 
rithmic singularities found in the S = 1 SU(2) case and 
extensions, although the physics of these anomalies can 
be expected to be similar. This discussion extends nat- 
urally to the underscreened lattice, and deserves clearly 
further theoretical studies. We think that a next step 
could be taken by solving the S = 1 Kondo lattice model 
in a combination of DMFT [2| and NRG techniques, and 
would help both to confirm the generic insulating behav- 
ior found here, and investigate in more detail the physical 
properties that are associated to it. 



2. Physical observation of the underscreened Kondo effect 

We finally discuss the possible experimental realization 
of underscreened Kondo physics. Concerning the single 
impurity case, the stability of S = 1 for quantum dots 
in even valleys of Coulomb blockade was discussed both 
theoretically and experimentally (for a review, see H^), 
although no emphasis was put onto the singular behav- 
ior that we have discussed previously. Another way to 
observe Kondo underscreening would be to study large S 
atoms deposited on a metallic surface and scanned by an 
STM tip. 

Finding an /-electron material related to the under- 
screened Kondo lattice seems to be a more daunting task. 
Indeed, not only the intersite magnetic correlations would 
need to be small in such a system, but a mechanism to 
form large spin (or spin-orbital degeneracy) of the /-level 
should be invoked. Typically, crystal fields leave unfortu- 
nately an effective Kramers doublet below temperature 
of several dozens of Kelvins. However, we think that the 
present work is interesting since it highlights the peculiar 
physical properties of underscreened lattices, and may be 
useful in strengthening an eventual experimental obser- 
vation. In fact, very recent transport measurements of a 
new cerium based compound, CeCuAs2, revealed weakly 
insulating behavior of the resistivity, with a power law 
divergence at low temperature |30| . While we are careful 
in connecting hastily this observation to our theoretical 
findinds, we hope that these remarks would stimulate fur- 



ther studies of underscreening effects in Kondo systems. 

V. CONCLUSION 

The underscreened Kondo model was investigated in 
this paper by means of a specially developed large N 
technique. Although the strong coupling fixed point in 
which itinerant electrons are tightly bound to the un- 
compensated spin is stable in the renormalization group 
sense, we have shown that non Fermi liquid behavior ap- 
pears in the form of anomalous power laws in the physical 
observables. The universal exponent was computed an- 
alytically and checked over the numerical solution of in- 
tegral equations, which show an interesting duality with 
respect to the previous theory of overscreening in mul- 
tichannel models. The extension to a finite dimensional 
lattice of underscreened magnetic moments was also con- 
sidered, and a pseudogap (weakly insulating) regime was 
discovered, with similar power laws governing the dy- 
namics of the electrons. We have also tried to point out 
several directions for future research, both concerning the 
theoretical aspects of underscreened models and their ex- 
perimental realization, for which CeCuAs2 might be a 
candidate. 
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APPENDIX A 

We present here the derivation of the non-trivial expo- 
nent governing the low frequency behavior of the physical 
quantities, following [l8ll24| . Let us assume the ansatz: 

TmG c {u) = A c \uj\- a " (Al) 
lmG B (uj) = A B \oj\~ aB Sgn{cj) (A2) 

at zero temperature. Using the spectral decomposition: 

G c (t) = / — e-^lmGciuj) (A3) 
Jo n 
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and the expression for the self-energies (|17I18|) with 
Gfo(r) = — |Sgn(r), we find: 

Ec(T) = _^£ r(1 _ afl) pi (A4) 
E.M = -^r(i-a =)R L_ 

Going back to frequency, we have simply: 



XmE c (w) = -^y?-\uj\ 



(A5) 



(A6) 



ImE fl (w) = -^M^SgnH (A7) 

To determine the real part in the previous self-energies, 
we use an analyticity argument, which gives for complex 



frequency z: 



i(l+a B )7r/2 



2 sin[(l + a s )7r/2]' 



e ia c 7r/2 

Ejj(z) = r-rk| 

K ' 2 sin[a c 7r/2]' 1 



(A8) 
(A9) 



and similarly for G c (z) and Gb(z). Finally, from Dyson's 
equation l|15 [l -l|lfi [l . we have G c {z) ~ — 1/Y, c (z) and 
Gb{z) ~ -1/Eb(z), providing relations between ampli- 
tudes A c , Ab and exponents a c , a^. After some manip- 
ulations, we find: 



2 1 
— arctan 

7T ^7 



(A10) 
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